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摘 要 : 区 域 土壤 盐分 时 空 变异 性 大 ,采用 经 典 统 计 和 地 统计 方法 
致 的 土壤 盐分 的 时 空 变化 趋势 。 本 文 以 内 蒙古 河套 灌区 隆 胜 研究 区 
分 数据 为 基础 ,利用 时 空地 统计 方法 分 析 区 域 土壤 盐分 时 空 变化 特征 ,比较 时 空 克 里 金 较 传统 空 
精度 提升 效果 ,并 验证 时 空地 统计 方法 在 监测 点 减少 50% 情 
围 是 0.43~1.14 ,为 中 强 变异 
积 盐 .农田 土壤 盐 


研究 区 土壤 盐分 空间 变异 系数 的 变化 范 
0.6-1.8 m 土壤 剖面 生育 期 脱盐 EAE 


结果 一 致 性 较 高 ,0~0.6 m 和 0.6~1.2 m 土壤 盐分 面 
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准确 判断 取样 时 间 不 规则 2s 
[X 68 个 监测 点 0~1.8m 土壤 剖面 4582 个 土壤 盐 


B 工 程 科 学 国家 重点 
呼和浩特 “010018 ) 


间 位 置 不 一 


青 况 下 预测 区 域 盐分 时 空 动态 的 能 力 。 
0-0.6 m 根 系 层 生 育 期 积 盐 、 非 生育 共 
mon (2) 和 度量 模型 
分 时 空 经 验 半 方 差 ,各 层 土壤 盐分 预测 值 和 观测 值 间 的 均 方 根 误差 RMSE 均 小手 0.21 dS m^ , 较 传 统 空 
RMSE 小 0.02~0.09 dS ^m^. (3) 采用 该 方法 在 减少 530% 监 疯 
PATA BAR 


上 点 情况 下 确定 的 土壤 盐 渍 化 分 布 与 所 有 取样 点 确定 的 
目 对 误差 MRE 分 别 为 -13.20% 和 -8.35% , RMSE 为 466.67 


s 间 克 里 金 插值 的 
结果 表明 :(1) 该 
脱盐 ， 
能 较 好 拟 合 盐 

s [B] yz FLA RJ 


hm 和 494.43 hm ,决定 系数 尼 为 0.79 和 0.72。 时 空 克 里 金 同 时 利用 了 土壤 盐分 时 间 和 空间 上 的 更 多 信息 ,实现 了 
稀 玻 盐分 监测 点 数据 集 土 壤 盐分 时 空 动 态 的 精确 估计 ,可 极 大 提高 区 域 土壤 盐分 时 空格 局 监测 的 效率 。 
关键 词 : 土壤 盐分 ; 时空 克 里 金 ; 河套 灌区 ; 时 空格 局 ; 和 度量 模型 


内 蒙古 河套 灌区 是 中 国 最 大 的 一 首 制 自流 灌 
盐碱地 占 全 部 灌溉 耕地 (57x10' hm?) 的 50% LA 
E ^, EHE Le Fe Ob. D TRE BEAR Mr Jr TT AY E 
要 障碍 ”。 精 确 评估 和 预测 土壤 盐分 时 空 分 布 , 实 
现 农 田 盐分 的 精准 管理 和 调控 对 提高 农业 生产 力 
意义 重大 。 然 而 ,河套 灌区 土壤 及 地 下 水 状况 复 
杂 ,加 之 强烈 的 人 类 活动 的 影响 ,土壤 盐分 时 空 变 
异性 强 ,实现 其 时 空 动态 预测 十 分 困难 。 学 者 们 一 
直 致 力 于 用 经 典 统计 学 和 地 统计 学 来 研究 土壤 盐 
分 的 变异 ,并 取得 了 一 些 成 果 ”。 在 内 蒙古 河套 灌 
区 , 徐 英 写 基 于 非 他 和 带 土壤 垂 向 水 分 和 盐分 的 时 
空 变异 性 ,用 条 件 模拟 理论 再 现 了 沙 壕 试 验 区 季节 
性 冻 融 十 水 盐 的 波动 过 程 ; 陈 亚 新 等 “研究 了 稳健 
NN en a A 
指示 克 立 格 与 神经 网 络 技术 融合 ,用 于 土壤 水 盐 
间 变 异性 评价 ; 史 海 滨 等 * 用 地 统计 学 方法 对 比分 
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分 布 特征 。 土 壤 盐分 是 具有 时 空 
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分 时 空 动 态 方面 ,Douaik 等 2 研究 了 匈牙利 Great 
Hungarian 平 原 一 块 25 hm? HE 0-0.4 m 土 层 的 土壤 
盐分 时 空 变化 ;Gasch 等 利用 和 度量 模型 建立 了 
美国 华盛顿 普尔 曼 附 近 一 块 37 hm 农田 3 维 土壤 电 

导 率 动态 模型 。 以 上 2 个 土壤 盐分 时 空 克 里 金 的 研 
究 均 是 田间 尺度 ,区 域 尺度 盐分 变异 性 增强 ,时 空 
克 里 金 方 法 在 区 域 土壤 盐分 时 空 建 模 与 预测 中 的 
应 用 效果 有 待 深入 人 研究 。 

本 文 以 内 蒙古 河套 灌区 隆 胜 研究 区 为 典型 研 
究 区 ,基于 2017 一 2018 年 研究 区 内 uad 
1.8 m 土壤 剖面 4582 个 土壤 盐分 数据 ,利用 时 空 
统计 方法 分 析 区 域 土 壤 盐分 时 空 变化 特征 ， Wee 

空 克 里 金 较 传统 空间 克 里 金 插 值 的 精度 提升 效果 ， 
并 验证 时 空地 统计 方法 在 监测 点 大 幅 减少 情况 下 
获取 区 域 盐 分 时 空 动态 的 能 力 ,人 研究 成 果 可 为 土壤 
盐 碱 化 时 空 动态 建 模 与 预测 提供 理论 和 科学 依据 。 


1 材料 与 方法 


1.1 研究 区 概况 

河套 灌区 深 处 内 陆 , 属 于 中 温带 干旱 半 干 旱 大 
陆 性 气候 ,冬季 严寒 少 雪 ,夏季 高 温 干 旱 。 隆 胜 研 
究 区 是 本 人 研究 的 主要 监测 区 ,位 于 内 蒙古 河套 灌区 
水 济 灌 域 ,西边 为 永 济 干渠 , 北 为 永 刚 分 干 沟 , 东 至 
东 济 文 渠 , 南 临 永 刚 分 干渠 ,西南 东北 长 约 15.5 km, 
西北 东南 宽 约 8.0 km, 总 土地 面积 8219.75 hn (图 
1)。 隆 胜 研究 区 每 年 11 月 中 旬 土 壤 开始 封冻 ,次 年 
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内 蒙古 河套 灌区 


图 1 隆 胜 研究 区 及 盐分 取样 点 位 置 图 
Fig. 1 Longsheng study area and soil salinity sampling 


locations 


一 层 的 土壤 水 盐 进 行 了 4 次 观测 ,分 别 标记 为 
Y1705 、Y1709 、Y1805 和 Y1809。 春 季 部 分 观测 点 地 


5 月 上 旬 融 通 。 邻 近 的 临河 气象 站 年 均 降 雨量 
148.8 mm, Æ 35 28 /z Œ (20 cm 28  IlIL) 2327 mm。 
2017 年 .2018 年 降雨 量 分 别 为 100.5 mm 、176.2 mm, 
生育 期 降雨 量 分 别 为 53.1 mm ,156.6 mm. 

研究 区 地 质 构 造 为 河 湖 相 交替 沉积 形成 的 湖 
相 和 河 相 沉积 层 ,0~1.0 m 土壤 质地 以 粉 砂 质 壤土 、 
粉 土 及 砂 质 壤土 为 主 , 水 平方 向 土壤 质地 分 布 极 不 
均匀 ,空间 变异 强烈 。0~2.5 m 垂 直方 向 上 土 体 构 型 
较为 复杂 ,多 有 黏土 AND Ie DB A EE. 
E d RES. 
12 土壤 盐分 采样 及 测试 

2017—2018 年 在 研究 区 内 均匀 布置 了 68 个 农 
田 土 壤 水 盐 观 测 点 (图 1) ,每 个 观测 点 距离 900 m Æ 
fio 分別 隆 2017 年 5 月 初 、2017 年 9 月 末 、2018 年 5 
月 初 、2018 年 9 月 末 対 0-1.8 m 深 度 范围 内 每 0.2 m 


下 水 埋 深 浅 , 取 到 地 下 水 面 处 停止 取样 。 每 个 观测 
点 用 土 钻 取 2 个 重复 ( 孔 ), 共 获取 4582 个 盐分 有 效 
样本 。 土 壤 盐 分 用 电导 率 仪 (上 海 雷 磁 DDSJ-308F) 
测试 水 土 比 为 5:1 的 土壤 浸 提 液 的 电导 率 。 文 中 所 
用 盐分 值 为 2 个 重复 的 均值 。 

图 1 中 的 32 个 长 期 监测 点 是 根据 改进 时 间 稳 
定性 方法 确定 的 ,将 其 用 以 验证 时 空地 统计 方法 在 
监测 点 减少 情况 下 获取 区 域 盐分 时 空 动态 的 能 
具体 确定 方法 见 文献 [24]。 

1.3 时 空地 统计 

时 空 克 里 金 方法 可 用 来 评估 土壤 盐分 的 时 空 
动态 。 通 常 普通 时 空 克 里 金 有 以 下 4 个 步 又 中 
(1) 时 空 经 验 半 方 差 函 数 的 计算 ; (2) 用 更 论 时 空 
半 方 差 拟 合 经 验 半 方差 ;(3) 用 时 空 克 里 金 做 时 空 
预测 ; (4) 用 留 一 交叉 验证 方法 评估 预测 结果 。 各 
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步骤 的 详细 计算 方法 如 下 : 
用 时 空 经 验 半 方差 水 数 y(hs, hn) 描述 土壤 参数 
时 空 分 布 的 相关 属性 1. 


N(hs,h;) 


Vhsh)= zrg gy Ly octhst hn) Zu] 
(1) 


SOP: hs JI PEAS ei a [IR] AT SEE ES BC Jes EP ES (m) ; hr 
JJ AVES js SE RT Jr ER E (d); Mhs, h) e PY 28 [8] ER 
(hs,h) 时 所 有 观测 样 点 的 对 数 ;Zs;,t) 和 2Z(si+hs,ti+h7) 
为 土壤 盐分 在 时 空位 置 (5;,#) 和 (siths,tithj) 的 值 。 

用 和 度量 模型 拟 合 时 空 经 验 半 方差 y(hs,h)"， 
该 模型 将 y(hs ,有 7) 分 为 空间 方差 ys(hs)、 时 间 方 差 yx(h7) 
和 联合 方差 yw(hs) ,三 者 关系 表示 为 : 

Vhs hy) my (hs) + Yrlhr)+Ysr(hsr) (2) 

式 中 eys(hs) yo) Fllyso(Asr) EY A BTA 388 A FH E 
间 理 论 半 方 差 函 数 模型 如 球状 模型 .高 斯 模型 、 指 
数 模型 来 表示 。 在 联合 方差 ysr(hsj) 中 第 引入 几何 各 
项 异性 比 a 来 计算 时 空 浪 后 距离 hs: 
hop = Jh} (ahy (3) 

该 研究 区 土壤 盐分 的 3 个 独立 方差 ys(pa,yrj) 
和 7ysr(ps 均 用 球状 模型 来 拟 合 。 球 状 模型 的 定义 如 
yp 


0 ; h,=0 
3h, h . 
y,(h) = ard E-E) O«h,xa, (i=S/TIST) 
C, * C; ; h,>a, 
(4) 
式 中 :Co 为 块 金 值 ; C 为 拱 高 ;Co+C 为 基 台 值 ;oa 为 相 


AH HE (ER LR s] .时间 和 时 空 变量 。 人 参照 
空间 相关 度 的 概念 ,用 各 部 分 独立 半 方 差 的 块 金 值 
之 和 与 基 台 值 之 和 的 比值 ($7D) 来 评价 时 空 相 关 
JE. STD«0.25 表示 时 空 相关 性 很 强烈 ,STD>0.75 
表示 时 空 相关 性 较 弱 ,其 余 表示 中 等 强度 时 空 相关 


性 


克 里 金 (Kriging) 插 值 是 对 区 域 变 量 取 值 进行 最 
优 线性 无 俩 估计 的 方法 。 时 空 克 里 金 插值 类 似 于 
空间 克 里 金 ,一 个 栅 格 中 的 土壤 盐分 值 ZrGo, 妇 由 邻 
近 的 取样 点 盐分 值 (s;,) 进 行 估算 


Z*(soto) = AZ の) (5) 


式 中 :为 已 知 观测 点 对 佑 值 点 的 加 权 系 数 。 根 据 
无 侦 佑 计 和 方差 最 小 两 项 要 求 可 确定 加 权 系 数 ”， 


21] 
o 


Zee en] FH A PIERR BER ASAE SEC 


[Y hsh) Yilhsshr) o Yohsshy) 1A 
Yalhshr) Yhsshr) tt Yalhsshr) 11A: 


y (hs. hr) 11^, 
OK 


Yahsh,) Yo 和 SP) … 
1 1 . 1 


和 (6) 
Violhs,hy) 


Ya fs hi) 
Y oh s. hr) 
| 1 
SU : yilhs, ho) J Z DRE Jes ER ES 29 hs RET [RETE Js EB 
H hi S Asi ti) (s t) EE 28 2E 77 288 se Bo] ME 
处 理 时 的 拉 格 朗 日 乘 数 ;为 搜索 域内 的 观测 数 。 

留 一 交叉 验证 用 来 评估 时 空 克 里 金 插 值 效 
果 。 交 叉 验证 时 ,验证 点 的 观测 值 被 移出 数据 集 ， 
用 剩余 的 数据 来 预测 该 时 空 点 的 土壤 盐分 值 呈 。 
然后 ,比较 土壤 盐分 观测 值 和 预测 值 的 拟 合 优 度 来 
评估 时空 克 里 金 模 型 。 拟 合 优 度 评价 指标 用 平均 
相对 误差 MRE、 均 方 根 误差 RMSE 和 決 定 系 数 だ , 
具体 计算 方法 见 文献 [24]。 

本 研究 基于 Matlab 环境 ,采用 遗传 算法 估计 变 
异 函 数 参 数 拉 ,并 完成 试验 半 方 差 散 点 图 、 理 论 模 
型 拟 合 曲面 图 的 绘制 ;时 空 插值 结果 分 布 图 在 Arc- 
GIS 10.3 中 完成 。 土 壤 盐 分 的 经 典 统计 分 析 利 用 
Excel 2019 软件 完成 。 


2 结果 与 分 析 


2.1 土壤 盐分 季节 性 变化 特征 

由 表 1 可 知 , 不 同 土 层 土壤 盐分 的 空间 均值 在 
0.25-0.35 dS+m! Z fH] 48 46 ,平均 处 于 非 盐 碱 化 水 
平 "。 土 壤 盐 分 的 峰 度 变化 范围 是 0.71~18.52, 偏 
度 变 化 范围 为 1.05~4.18 ,数据 右 侧 拖 尾 严重 ,所 有 
取样 时 间 各 层 土壤 盐分 均 是 非 正 态 分 布 的 。 土 壤 
盐分 空间 变异 系数 的 变化 范围 是 0.43~1.14 ,为 中 强 
变异 。0~0.4 m 越 接近 表层 ,土壤 盐分 受降 雨 灌溉 、 
蒸发 等 上 边界 作用 越 强烈 ,盐分 随时 间 波 动 越剧 
烈 。0.4-0.6 mm 土壤 盐分 均值 不 同 观测 期 变化 范围 
为 0.31~0.32 dS -m ,变化 幅度 不 超过 5% ,是 盐分 最 
为 稳定 的 土屋 。2017 年 .2018 年 生育 初期 及 生育 末 
期 0~0.6 m 根 系 层 土 壤 盐 分 分 别 为 0.32 dS +m” ,0.34 
dS*m^,0.29 d$: m^,0.30 dS.m"', 根 系 层 土壤 生育 期 
处 于 积 盐 状态 。 以 2017 年 生育 期 末 和 2018 年 生育 
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表 1 4 次 取样 时 间 土 壤 盐 分 统计 特征 
Tab.1 Statistical characteristics of soil salinity at four sampling times 


均值 中 位 数 标准 差 变异 峰 度 備 度 最小 値 RAE 


土 层 /m 时 间 观测 数 


/(dS:m') /(d$-m?)  /(d$*m?) 系数 /(d$.m) /(dS*m*) 
0-0.2 Y1705 65 0.34 0.23 0.39 1.14 18.52 4.18 0.13 2.34 
Y1709 66 0.35 0.23 0.28 0.78 Sil 2.12 0.1 155 
Y1805 67 0.28 0.23 0.13 0.45 3.34 1.69 0.14 0.77 
Y1809 66 0.28 0.2 0.21 0.77 11.34 3.18 0.08 1.29 
0.2~0.4 Y1705 65 0.3 0.25 0.18 0.6 3.84 1.87 0.12 0.99 
Y1709 66 0.33 0.24 0.24 0.74 3.39 1.97 0.12 1.19 
Y1805 67 0.29 0.23 0.15 0.52 1.39 1.42 0.13 0.76 
Y1809 66 0.3 0.23 0.21 0.71 9.18 2.87 0.08 1.23 
Y1705 65 0.31 0.27 0.21 0.66 10.29 2.74 0.11 1.35 
0.40.6 Y1709 66 0.32 0.26 0.2 0.63 1.99 1.57 0.07 0.96 
Y1805 67 0.32 0.24 0.2 0.63 2.58 1.64 0.12 1.02 
Y1809 66 0.32 0.26 0.2 0.63 3.14 1.77 0.1 1.03 
0.6~0.8 Y1705 65 0.33 0.26 0.2 0.6 4.3 1.88 0.1 1.18 
Y1709 66 0.29 0.26 0.15 0.53 1.27 1.21 0.07 0.74 
Y1805 67 0.33 0.26 0.2 0.61 2.22 1.62 0.12 1.03 
Y1809 66 0.31 0.24 0.2 0.64 3.73 1.87 0.09 1.06 
0.8~1.0 Y1705 65 0.33 0.27 0.19 0.57 5.64 2.06 0.1 1.18 
Y1709 66 0.28 0.27 0.13 0.47 1.94 1.17 0.07 0.73 
Y1805 67 0.33 0.27 0.18 0.54 2.17 1.53 0.12 0.92 
Y1809 66 0.29 0.22 0.17 0.61 3.3 1.8 0.09 0.89 
1.0-1.2 Y1705 64 0.31 0.26 0.16 0.52 3.7 1.76 0.06 0.9 
Y1709 66 0.28 0.25 0.13 0.46 0.71 1.1 0.09 0.63 
Y1805 67 0.33 0.27 0.17 0.53 2.26 1.53 0.11 0.9 
Y1809 66 0.28 0.23 0.18 0.62 8.86 2.63 0.09 113: 
1.2-1.4 Y1705 $1 0.29 0.24 0.15 0.51 5.79 1.98 0.12 0.92 
Y1709 66 0.26 0.24 0.12 0.44 1.83 1.24 0.09 0.68 
Y1805 65 0.31 0.25 0.15 0.48 0.78 1.26 0.11 0.72 
Y1809 66 0.27 0.23 0.15 0.56 4.2 1.85 0.1 0.91 
1.4-1.6 Y1705 45 0.28 0.23 0.15 0.54 4.49 1.82 0.12 0.88 
Y1709 66 0.25 0.24 0.11 0.43 0.89 1.05 0.1 0.55 
Y1805 64 0.29 0.24 0.13 0.46 2.4 1.53 0.14 0.78 
Y1809 66 0.26 0.22 0.13 0.51 6.12 2.17 0.08 0.86 
1.6-1.8 Y1705 29 0.32 0.27 0.19 0.61 5.14 2.04 0.12 1.01 
Y1709 66 0.26 0.24 0.12 0.48 5.72 1.76 0.07 0.81 
Y1805 58 0.27 0.22 0.12 0.45 3.16 1:52 0.11 0.77 
Y1809 66 0.25 0.23 0.12 0.47 2.99 1.64 0.08 0.65 


HIIR E ERRIA BABE TT PT, BATE AR AE 
用 使 根系 层 土壤 盐分 降低 了 14.7% 。 这 与 孙 贯 芳 
等 ?9 在 河套 灌区 秋 浇 灌溉 洗 盐 试验 的 结果 基本 一 
致 : 秋 浇 灌 黄 河水 180 mm 后 ,次 年 春播 前 0~1.0 m 
土壤 盐分 平均 下 降 10.86%~26.14%。 从 整个 剖面 来 
看 ,5 月 初 土壤 盐分 0.6~1.0 m 土 层 盐分 最 大 ,两 端 较 


小 , 秋 浇 及 冻 融 作用 使 盐分 在 根系 层 底 部 以 下 剖面 
聚集 ,生育 期 蒸发 蒸腾 作用 又 使 部 分 盐分 向 上 运 
动 ,导致 0.6m 以 下 各 层 土壤 盐分 生育 初期 明显 大 
于 生育 末期 。 忽 略 2018 年 科大 降雨 対 0-0.4 m E 
盐分 的 淋 洗 作用 ,整个 土壤 剖面 9 月 底 盐分 分 布 自 
上 到 下 依次 减 小 ,属于 典型 的 灌溉 蒸发 型 土壤 盐分 
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前面 。 2.3 土壤 盐分 时 空 分 布 的 交叉 验证 与 预测 


2.2 土壤 盐分 时 空 半 方差 

用 对 数 变换 后 的 0~1.8 m 土壤 盐分 数据 计算 经 
验 时 空 半 方差 ,结果 如 图 2。 所 有 半 方 差 在 变 程 范 
围 内 均 随 着 时 间 沾 后 距离 和 空间 滞后 距离 的 增 大 
而 增 大 ,说 明 土 壤 盐 分 的 相似 性 随 着 时 间 和 空间 距 
离 的 增 大 而 减 小 。 土 壤 盐分 的 时 空 半 方差 值 逐 渐 
达到 一 个 常数 值 ,对 数 转换 后 的 土壤 盐分 数据 在 空 
间 上 是 均匀 的 ,并 具有 时 间 稳 定性 ””"。 用 和 度量 
模型 拟 合 经 验 时 空 半 方差 ,并 对 所 有 的 时 空位 置 重 
复 进 行 交 叉 验 证 ,基于 交叉 验证 结果 选择 的 时 空 半 
方差 模型 如 表 2 所 示 。 通 常 很 难 区 分 空间 .时间 和 
联合 部 分 的 块 金 值 , 表 中 参数 CG, 是 3 部 分 块 金 值 
和 ,不 同 土 层 土壤 盐分 块 金 值 的 变化 范围 为 0.0028~ 
0.0186, 基 台 值 变化 范围 为 0.0371~0.0618。0.4~0.6 
m 土 层 的 时 空 相关 度 STD 等 于 0.32, 属 于 中 等 相关 ， 
其 余 的 时 空 相 关 度 变化 范围 为 0.06~0.18, 均 小 于 
0.25 ,时 空 相关 性 强烈 。 


(a) 0~0.2m 


“200 _ 
wa 


(b) 0.2-0.4 m 


土壤 盐分 时 空 克 里 金 交 叉 验 证 结果 如 图 3 所 
示 , 各 层 土壤 盐分 预测 值 和 观测 值 间 的 MRE 变化 范 
围 为 3.38%~5.59% 。0~0.2 m 土壤 盐分 的 RMSE 是 
0.21 dS- m” ,其 他 层 土壤 盐分 RMSE 均 小手 0.15 dS- 
m'。 各 层 土 壤 盐 分 预测 值 和 观测 值 间 的 决定 系数 
尼 在 0.50~0.66 间 变化 。 该 土壤 盐分 交叉 验证 的 尼 
结果 远 好 于 Gasch 等 报道 的 土壤 盐分 时 空 克 里 金 
交叉 验证 结果 0.18, 基 本 上 接近 于 Douaik 等 2 报道 
的 土壤 盐分 时 空 克 里 金 验证 结果 0.76。 由 图 3 可 
知 , 除 少数 土壤 盐分 值 较 大 而 预测 结果 偏 小 的 情况 
外 ,大 部 分 预测 的 土壤 盐分 与 实测 数据 吻合 良好 ， 
和 度量 模型 能 够 合理 表征 土壤 盐分 的 时 空 结 构 ,并 
取得 精确 的 预测 结果 。 

此 外 ,传统 空间 克 里 金 方法 计算 的 各 时 期 半 方 
差 模型 见 文献 [24]。 土 壤 盐分 的 空间 结构 受 自 然 
和 人 类 活动 的 影响 ,4 个 取样 时 间 有 较 大 的 差别 。 
以 0.4~0.6 m + Je X ffl , Y1705, Y1709, Y1805, 


图 2 各 层 经 验 ( 蓝 色 圆圈 ) 和 拟 合 的 土壤 盐分 时 空 半 方 差 曲 面 


Fig.2 Empirical (blue circle) and fitted spatio-temporal semivariogram surfaces of the soil salinity in each layer 
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表 2 土壤 盐分 时 空 半 方差 函数 模型 及 交叉 验证 结果 


Tab. 2 The semivariogram models and parameters of soil salinity 


PM o 空间 参数 时 间 参 数 联合 参数 
Cs as Cr ar Csr asr a 

0~0.2 0.0063 0.0240 957.00 0.0282 292.70 0.0032 235.01 0.1000 
0.2-0.4 0.0081 0.0143 1438.69 0.0201 287.37 0.0092 1217.14 0.6647 
0.4-0.6 0.0186 0.0152 1001.61 0.0073 289.24 0.0162 1868.67 0.1000 
0.6-0.8 0.0080 0.0249 1615.07 0.0147 227.55 0.0074 1003.57 0.1856 
0.8~1.0 0.0055 0.0259 1916.53 0.0188 209.64 0.0009 1003.37 0.4790 
1.0~1.2 0.0028 0.0257 1323.09 0.0158 220.30 0.0009 955.91 3.0303 
1.2~1.4 0.0069 0.0169 863.06 0.0061 147.70 0.0093 1930.13 0.2737 
1.4~1.6 0.0061 0.0187 1127.09 0.0063 175.33 0.0061 1879.02 1.1730 
1.6~1.8 0.0042 0.0010 1863.53 0.0132 293.88 0.0241 1827.54 0.0978 
0~0.6 0.0089 0.0280 1250.41 0.0046 246.15 0.0017 1861.18 0.0244 
0.6-1.2 0.0054 0.0290 1945.91 0.0090 244.73 0.0013 1259.28 0.0391 


Y1809,4 个 时 期 变 程 分 别 为 1610 m、1680 m 、1801 m 
和 6180 m, 蔡 取样 时 间 刚 好 是 Y1809, 将 会 错误 的 认 


局 。 该 研究 区 土壤 盐分 空间 分 布 格局 受 地 势 E 
质地 下 水 埋 深 及 矿 化 度 等 多 种 因素 的 影响 ,人 研究 
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为 该 区 域 0.4~0.6 m 土壤 盐分 的 变 程 是 6000 m Æ 
^ , 远 超 实际 变 程 1600 m, 用 Y1809 取得 的 结果 认 
识 区 域 土壤 盐分 空间 结构 特征 ,将 存在 严重 的 偏 
差 , 仅 由 1 次 取样 获取 土壤 盐分 的 空间 结构 是 不 可 
靠 的 。 因 此 , 需 多 次 监测 综合 分 析 确 定 区 域 土 壤 盐 
分 的 空间 结构 特征 。 时 空 克 里 金 能 够 同时 利用 不 
同 监测 时 间 的 盐分 信息 ,是 解决 这 一 问题 的 有 力 手 
段 。 图 4 展示 了 传统 空间 克 里 金 和 时 空 克 里 金 方法 
计算 的 各 时 期 各 土 层 间 土 壤 盐 分 交叉 验证 的 RMSE 
值 ,时 空 克 里 金 的 RMSE 较 传统 空间 克 里 金 小 0.02~ 
0.09 dS.m', 显 著 提高 了 土壤 盐分 预测 精度 。 
基于 半 方 差 函 数 分 析 , 用 时 空 克 里 金 插值 方法 
评估 了 不 同 取 样 时 间 各 层 土 壤 盐 分 的 空间 分 布 。 
土壤 盐分 空间 分 布 具有 时 间 稳 定性 ,不 同时 间 盐 分 
空间 分 布 格局 基本 一 人 致 池 。 图 5 展示 了 Y1805 各 层 
土壤 盐分 的 空间 分 布 。 同 一 时 间 不 同 层 间 土壤 盐 
分 分 布 具有 相似 性 ,各 层 土壤 盐分 高 低 值 分 布 位 置 
较为 接近 ,特别 是 相 邻 层 间 盐分 斑 块 位 置 较 为 一 
致 , 仅 斑 块 的 大 小 和 范围 有 所 差别 ,说 明 土 壤 盐分 
垂 向 分 布 较为 连续 ,各 层 土壤 盐分 含量 密切 相关 。 
0~0.6 了 各 层 土壤 盐分 高 低 相 间 , 呈 插 花 式 分 布 ,总 
体 表现 为 研究 区 东部 及 东北 部 盐分 含量 高 ,而 西北 
部 和 南部 盐分 含量 较 低 。0.6~1.8 m 各 层 土 壤 盐 分 
在 研究 区 中 部 和 北部 沿 东 北 -西南 方向 呈 条 人 带 状 分 
布 , 自 西向 东 表 现 为 低 -高 - 低 - 高 的 盐分 分 布 格 


区 东部 排水 沟 经 过 位 置 处 盐分 含量 反而 较 其 他 位 
置 明显 大 ,主要 是 由 于 该 区 域 地 势 低 , 地 下 水 埋 深 
浅 矿 化 度 高 , 且 土 壤 多 为 粉 土 和 黏土 汪 ,排水 沟 年 
久 失 修 未 能 有 效 发 挥 作用 所 致 。 
2.4 时 空地 统计 获取 稀 疏 监测 点 盐分 时 空 动态 的 
能 力 

0~0.6 m 根 系 戻 和 0.6-1.2 m 土 层 时 空 半 方差 模 
型 及 参数 见 表 2, 留 一 交叉 验证 结果 如 表 3。 长 期 监 
测 点 确定 的 0~0.6 m 和 0.6-1.2 m 土壤 盐分 预测 值 与 
观测 值 的 MRE 分 别 为 2.52% 和 2.54% , RMSE 为 
0.12 dS・m! 和 0.09 ds.m-, 尼 为 0.73 和 0.70。 长 期 
监测 点 确定 的 根系 层 土壤 盐分 空间 分 布 格局 (图 
6al~6d1) 和 所 有 取样 点 确定 的 空间 分 布 图 (图 6a2~ 
6d2) 基本 一 致 。 如 图 7a .图 7c 所 示 ,由 土壤 盐分 长 
期 监测 点 确定 的 不 同 盐 碱 化 水 平 的 面积 与 所 有 取 
样 点 确定 的 面积 吻合 较 好 ,两 者 间 的 MRE 为 
-13.20% ,RMSE 为 466.67 hm’, R X 0.79 ,说 明 由 土 
壤 盐 分 长 期 监测 点 确定 根系 层 土壤 盐分 的 空间 分 
布 是 可 行 的 。 同 样 地 ,由 土壤 盐分 长 期 监测 点 确定 
的 0.6-1.2 m 土壤 层 盐 分 分 布 (图 6e2~6h2) 和 所 有 取 
样 点 确定 的 盐分 分 布 图 (图 6el~6hD) 也 高 度 一 致 。 
两 者 确定 的 相同 盐 碱 化 水 平 的 面积 统计 结果 也 较 
为 接近 (图 7b .图 7d, 面 积 间 的 MRE 为 -8.35% , 
RMSE 为 494.43 hm , 尼 为 0.72, 说 明 土 壤 盐 分 长 期 
监测 点 可 以 用 来 确定 0.6~1.2 m 土壤 层 盐分 的 空间 
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用 于 时 空 克 里 金 交叉 验证 的 序列 号 


图 3 土壤 盐分 时 空 克 里 金 交叉 验证 结果 


Fig.3 Soilsalinity cross-validation results for the spatio-temporal Kriging 
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图 4 传统 空间 克 里 金 和 时 空 克 里 金 土壤 盐分 交叉 验证 结 
果 比 较 
Fig. 4 The cross-validation results of traditional spatial 


Kriging and spatio-temporal Kriging 


分 布 。 XE— AG LOBE 2S vu HUS Sj TOR v HUE TROC 
iri UI e ERP NEA JEU GE 77 ,两 者 土壤 盐分 长 
期 监测 点 确定 的 不 同 盐 碱 化 水 平 的 面积 与 所 有 取 
样 点 面积 间 的 评价 指标 结果 见 表 4, 时 空 克 里 金 的 
结果 明显 优 于 普通 克 里 金 ,说 明 在 数据 监测 点 较为 
稀 玻 的 情况 下 , 相 较 于 普通 克 里 金 ,时 空 克 里 金 也 
提高 了 稀 玻 监测 点 土壤 盐分 时 空 动 态 的 预测 精度 。 
综 上 所 述 , 由 土壤 盐分 长 期 监测 点 确定 的 土壤 


(a) 0-0.2 m ; (b) 0.2~0.4 m ; (c) 0.4~0.6 m ; (d) 0.6-0.8 m 


盐分 空间 分 布 状况 与 所 有 取样 点 较为 一 致 ,利用 时 
空 克 里 金 在 土壤 盐分 监测 点 数据 减少 为 原来 一 半 
情况 下 仍 可 以 精确 估计 土壤 盐分 时 空 动态 ,以 较 少 
的 采样 数据 实现 了 土壤 盐分 时 空 变化 特征 的 精确 
评 佑 , 极 大 提高 区 域 土壤 盐分 监测 的 效率 。 


3 结论 


本 文 以 内 蒙古 河套 灌区 隆 胜 研究 区 为 典型 研 
究 区 ,基于 2017 一 2018 年 研究 区 内 68 个 监测 点 0~ 
1.8 m 土壤 盐分 数据 ,分 析 了 土壤 盐分 的 季节 性 变化 
特征 ,利用 时 空地 统计 方法 研究 了 区 域 土壤 盐分 时 
空 变化 特征 及 时 空 克 里 金 的 插值 精度 ,并 进一步 验 
证 时 空地 统计 方法 在 监测 点 不 足 原来 一 半 情 况 下 
获取 区 域 盐分 时 空 动 态 的 能 力 。 主 要 结论 有 : 

(1) 该 研究 区 土壤 盐分 有 明显 的 季节 性 规律 ， 
0~0.6 m 根 系 层 生 育 期 积 盐 , 非 生育 期 脱盐 ;0.6~1.8 
m 土壤 剖面 生育 期 脱盐 , 非 生育 期 积 盐 。 

(2) 和 度量 模型 能 较 好 拟 合 盐分 时 空 经 验 半 方 
差 , 时 空 友 里 金 对 土壤 盐分 交叉 验证 的 RMSE 较 传 


(e) 0.8~1.0m 


(i) 1.6-1.8 m EC/(dS-m) 


mz 


(f) 1.0~1.2 m j (g) 1.2~1.4m 3 T A 


图 5 0-1.8 m 各 层 土壤 盐分 空间 分 布 (Y1805) 
Fig. 5 Soil salinity maps of different soil layers within the depth 0-1.8 m in Y1805 


表 3 由 所 有 盐分 监测 点 (ASML) 和 长 期 监测 点 (LSML) 确 定 的 土壤 盐分 时 空 克 里 金 验证 结果 


Tab.3 Cross-validation results for 0-0.6 m and 0.6-1.2 m layer soil salinity with all soil salinity monitoring locations 


(ASML) and long-term soil salinity monitoring locations (LSML) 


ASML LSML 
土 层 /m 
MRE/% RMSE/(dS-m"') R MRE/% RMSE/(dS* m) 
0-0.6 3.60 0.13 0.64 2.52 0.12 0.73 
0.6~1.2 3.11 0.10 0.65 2.54 0.09 0.70 
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图 6 由 土壤 盐分 长 期 监测 点 LSML) 和 所 有 监测 点 (ASML) 确 定 的 0~0.6 m 根 系 层 和 0.6~1.2 mm 土壤 层 盐分 空间 分 布 比较 
Fig.6 The spatial distribution of root zone (within the depth of 0-0.6 m) and 0.6-1.2 m layer soil salinity determined by all soil 


salinity monitoring locations (ASML) and long-term soil salinity monitoring locations (LSML) at the 4 sampling times 


—o— ASML -—4—LSML 


(a) 0~0.6 m 
MRE- -13.20% RMSE-466.47 hm? R?=0.79 


3600. (b 0.61.2 m 


IMRE= -8.35% RMSE-494.43 hm? 
30004  g-0.72 
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图 7 由 土壤 盐分 长 期 监测 点 LSML) 和 所 有 取样 点 (ASMI) 确 定 的 不 同 土壤 盐 碱 化 面积 的 比较 


Fig.7 The area of different soil salinity determined by all soil salinity monitoring locations (ASML) and long-term soil salinity 


monitoring locations (LSML) 
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克 里 金 获取 稀 玻 监测 点 土壤 盐分 时 空 
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动态 能 力 的 比较 


Tab.4 Comparison of ability of spatio-temporal Kriging and ordinary Kriging to obtain soil salinity dynamics under sparse 


monitoring locations 


时 空 克 里 金 普通 克 里 金 
土 层 /m 
MRE/% RMSE/hm* R MRE/% RMSE/hn R 
0~0.6 -13.20 466.47 0.79 26.14 412.99 0.78 
0.6-1.2 -8.35 494.43 0.72 120.09 731.13 0.48 
BRAS 回 克 里 金 小 0.02~0.09 dS・ m’, 时 空 s 克 里 金 同 时 Ying, Wei Zhanmin, et al. Approximation method for spatial vario- 
利用 了 土壌 +h 分 时 间 和 空 间 E 的 更 を 信息 dy 3t xs 通 gram of water and salinity in soil based on robust statistics[J]. Jour- 
Hir FI AU 5 
EXE dz 、 Ir e" nal of Hydraulic Engineering, 2004, 35(9): 46-51. | 
克 里 金 显著 提高 土壤 盐分 的 预测 精度 。 Mcd . 
- " ES Be ENDL 、 [7] 刘 全 明 . 非 参数 统计 与 人 工 智能 技术 在 水 土 空间 变异 中 的 应 用 
(3) 同一 时 间 不 同 层 间 土壤 盐分 分 布 具有 相似 ——— 
研究 [D]. 呼和浩特 : 内 蒙古 农业 大 学 , 2009. [Liu Quanming. Ap- 
J AY EL ANT e mM RS £ 
性 ， 相 邻 士 层 间 盐 分 5 二 块 位 置 较为 x, 仅 斑 块 的 plication of Non- parametric Statistics and Artificial Intelligence 
大 小 和 范 iu 围 有 所 差别 。 0~0.6m 各 层 土壤 盐分 高 低 Technology in Spatial Variability of Water and Soil[D]. Hohhot: In- 
相间 , 呈 插 花 式 分 布 ,0.6~1.8 mm 各 层 土壤 盐分 在 研 ner Mongolia Agricultural University, 2009. ] 
究 区 中 部 和 北部 沿 东北 -西南 方向 呈 条 带 状 分 布 ， [8] SEHRE, 吴迪 , 同 建 文 , 等 . 盐 渍 化 灌区 节 水 改造 后 土壤 盐分 时 
空 变化 规律 研究 由 . 农业 机 械 学 报 , 2020, 51(2): 318-331. [Shi 
自 西 向 东 表现 为 低 - 高 - 低 - 高 的 盐分 分 布 格局 。 a 
し Haibin, Wu Di, Yan Jianwen, et al. Spatial-temporal variation of 
(4) EAS ERRETA ULT ict — 91 SIR th wn l MCN 
soil salinity after water saving transformation in salinized Irriga- 
人 hk ZA pA A 
测 点 土壤 盐 分 数据 ? 2h H 時 人 we 里 金 2 计算 的 中 Heth tion District|J]. Transactions of the Chinese Society for Agricultur- 
分 空间 分 布 状况 与 所 有 取样 点 结果 较为 一 至, 实现 al Machinery, 2020, 51(2): 318-331. ] 
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Spatio-temporal patterns of soil salinity in Hetao Irrigation District based on 
spatio-temporal Kriging 


SUN Guanfang, GAO Zhaoliang, ZHU Yan, YANG Jinzhong, QU Zhongyi 
(1. Institute of Soil and Water Conservation, Northwest A & F University, Yangling 712100, Shaanxi, China; 
2. State Key Laboratory of Water Resources and Hydropower Engineering Science, Wuhan University, Wuhan 
430072, Hubei, China; 3. College of Water Conservancy and Civil Engineering, Inner Mongolia Agricultural 


University, Hohhot 010018, Inner Mongolia, China) 


Abstract: For maintaining crop yield in salt-affected dry agricultural settings, monitoring and analyzing spatio- 
temporal dynamics of soil salinity over broad areas is crucial yet challenging due to its high variability. The most 
popular techniques for evaluating spatial distribution patterns and temporal trends are classical statistical analysis 
and traditional geostatistical analysis, but they are not suitable for accurately capturing spatio-temporal trends of 
soil salinity due to irregular sampling time and inconsistent spatial position during sampling time. Spatio- 
temporal Kriging is an extension of spatial geostatistics to space-time geostatistics and may overcome this 
problem effectively because its model covariance/variance is a function of both space and time. However, its 
application in spatio-temporal modeling and prediction of regional soil salinity is still unclear. Based on 4582 soil 
salinity data of 0— 1.8 m soil profiles from 68 monitoring locations in the Longsheng study area of Hetao 
Irrigation District, Inner Mongolia, spatio- temporal variation characteristics of regional soil salinity using a 
spatio- temporal geostatistical method, and spatio-temporal Kriging interpolation accuracy was compared with 
traditional spatial Kriging interpolation. Furthermore, the ability of spatio- temporal Kriging to obtain regional 
soil salinity dynamics was verified using less than half of the original monitoring locations. The results showed 
that the spatial variation coefficient of soil salt in the study area ranged from 0.43 to 1.14, which was categorized 
as medium to strong variability. Regional averaged soil salinity dynamics had obvious seasonal variation 
characteristics, and the root zone (0—0.6 m) soil salinity accumulated in the crop growing season and desalted in 
the fallow season, while the deep soil salinity (0.6-1.8 m) was the opposite. The sum- metric model can fit the 
temporal and spatial experience semi- variance of soil salinity well, and the root mean square error (RMSE) 
between the predicted value and observed value of soil salinity in each layer was less than 0.21 dS・m_ , which 
was 0.02-0.09 dS・m' less than that of traditional spatial Kriging. The areas of different soil salinity determined 
by 32 sparse monitoring locations were in good agreement with those determined by all sampling sites, and the 
mean relative error between areas of different soil salinity for 0-0.6 m and 0.6-1.2 m were —13.20% and —8.35%, 
respectively. Similarly, the respective RMSE were 466.67 hm’? and 494.43 hm’ and the determination coefficient 
(R^) were 0.79 and 0.72, indicating that spatial distribution of soil salinity obtained by sparse monitoring locations 
is consistent with the results of all sampling locations. Spatio- temporal Kriging significantly improves the 
prediction accuracy of soil salinity compared with ordinary Kriging, since it uses more information on soil 
salinity in time and space. The accurate estimation of spatio-temporal dynamics of soil salinity in the data set of 
sparse monitoring points was realized, which can greatly improve the monitoring efficiency of the spatio- 
temporal pattern of soil salinity in the region. 

Keywords: soil salinity; spatio-temporal Kriging; Hetao Irrigation District; spatio-temporal pattern; sum-met- 


ric model 


